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در مواجه با خطر سیل و یا خسارات ناشی از خشکسالی, برآورد میزان بارش و الگوی تغییرات مکانی آن 
در یک منطقه گسترده. یکی از چالش‌های مهم در علوم هواشناسی. کشاورزی و هیدرولوژی است. 
اندازه‌گیری محلی بارندگی در مناطق دور افتاده به دلیل هزینه زیاد و محدودیت‌های عملیاتی دشوار است. 
بدین علت در تحقیق حاضر به‌منظور تعیین الگوی مکانی-زمانی بارش و امکان تلفیق داده‌ها؛ سه نوع 
مختلف از تولیدات بارندگی شامل داده‌های ماهواره‌ای (118]135342), داده‌های حاصل از مدل پیش‌بینی 
عددی جوی (۷115) و اندازه گیری‌های زمینی (نقشه‌های حاصل از روش زمین‌آمار (1600))؛ مورد 
مطالعه قرار گرفتند. این مطالعه در بازه زمانی سال‌های ۲۰۰۰ تا ۲۰۱۰ میلادی و برای منطقه شمال شرق 
ایران به صورت ماهانه فصلی و سالانه انجام شد. داده‌ها با استفاده از شاخص اعتبارسنجی ۷155 و 
الگوریتم تشابه با یکدیگر مقایسه شدند. نتایج نشان دادند که یکی از ضعف‌های روش زمین‌آمار نبودن 
اطلاعات کافی در ارتفاعات بالای (۱۵۰۰) متر منطقه است. همچنین دقت تصاویر ماهواره‌ای در فصل‌های 
گرم بیشتر بود؛ بطوریکه در ماه آگوست مقدار ۱/۸۷ < 15۳8 به دست آمد. در فصل زمستان (ماه ژانویه) 
بیشترین مقدار ۱۶/۰۲ < 158 حاصل شد که این امر عملکرد ضعیف تولیدات ماهواره‌ای ]112۷۷ در 
مناطق پوشیده از یخ را نشان می‌دهد. در اعتبارسنجی مدل ۷1۳15 بیشترین و کمترین مقدار 1318۳ به 
ترتیب 1/۶ و ۱/۰۵ به دست آمد. علاوه بر این مدل 1۷۲5 تا حدود زیادی در شبیه‌سازی مقادیر بارندگی 


سالانه بیش‌برآورد داشت. نتایج تحلیل‌های مکانی- زمانی الگوریتم تشابه نیز نشان دادند که عملکرد مدل 


۱ نویسنده مسئول : ۰۵۱۳۸۸۰۵۷۱۷ مه اتقصصع ۵ صتعقمصعه :اتمصرظ 


۱۳۹۶ جغرافیا و مخاطرات محیطی شماره بیست و سوم 


5 در مقیاس ماهانه و فصلی و تعیین مناطق بارندگی بهتر از تصاویر ماهواره‌ای 11811 بود. همچنین 
هر سه محصول الگوی مکانی بارندگی در مقیاس فصلی و سالانه را به‌حوبی نشان دادند. 


۱. مقدمه 


تخمین دقیق بارندگی در مناطقی با ویژگی‌های جغرافیایی پیچیده در زمینه اقلیم‌شناسیء هواشناسی کشاورزی و 
هیدرولوژی از اهمیت ویزه‌ای برخوردار است؛ اما وسعت نواحی و تغییرات شدید زمانی -مکانی بارش 
اندازه‌گیری‌ها را مشکل می‌سازد. عدم وجود داده‌های مکانی بارش از چالش‌های مهم در پیش‌بینی وقوع سیل یا 
خشکسالی و هشدار به‌موقع در اين موارد است. یکی از راهکارهای ارائه شده در این مورد. اندازه‌گیری بارندگی از 
فضا است. تکنیک سنجش ازدور تصاویر کامل‌تری از ساختار سیستم‌های بارندگی ارائه می‌کند و در چند دهه اخیر 
توسعه چشمگیری یافته است. از این میان ماهواره ۳۲181۷1 اولین تلاش بین‌المللی برای اندازه‌گیری قابل اعتماد 
بارش از فضا می‌باشد (اسمیت " و همکاران ۲۰۰۷: ۱۱۱). یکی دیگر از مجموعه داده‌هابی که در سال‌های اخیر در 
دسترس قرار گرفته است. خروجی مدل‌های پیش‌بینی عددی است. مدل‌های مذکور با حل معادلات اولیه حاکم بر 
جو و زمین به محاسبه پارامترهای مختلف محیطی و هواشناسی مانند بارندگی در پیکره شبکه‌بندی شده کره زمین 
می‌پردازند (میج " و همکاران. ۲۰۰۹ :10۱۱). هر یک از مجموعه داده‌های ذکر شده نقاط ضعف و قوت خود را 
دارند. به‌عنوان‌مثال اندازه‌گیری‌های زمینی در عین اينکه تنها اندازه‌گیری مستقیم بارندگی هستند در مناطق کوهستانی 
پراکندگی مکانی اندکی دارند (اسچیل * و همکاران»۲۰۱۱: ۲04۹). باران‌سنج‌ها اغلب در مجاورت روستاها نصب 
می‌شوند بدین ترتیب در مناطق غیرمسکونی عملا نمونه اندازه‌گیری بارش وجود ندارد. از مزایای تولیدات 
سنجش زدور می‌توان به پوشش نواحی زیادی از سطح زمین اشاره کرد بااین‌حال تصاویر مربوط به اواسط دهه 
۰ که فقط با تکیه بر امواج مادون قرمز تهیه شده‌انده در تخمین بارندگی و مخصوصاً بارندگی‌های همرفتی و 
گرادیان شدید بارش دقت کافی را ندارند. محصولات جدید ماهواره‌ای از ترکیب امواج مادون‌قرمز و مایکروویو و 
داده‌های زمینی برای شناسایی بارش‌های شدید استفاده می‌کنند. اما هنوز تفاوت چشمگیری میان تولیدات مختلف 
ماهواره‌ای و استفاده از آن‌ها در مناطق مختلف وجود دارد (وارد* و همکاران. ۲۰۱۱: ۱۲۲۲). همچنین تولیدات 
ماهواره در تشخیص بارش مایع از جامد عملکرد ضعیفی دارند. گروه دیگری از داده‌ها در بررسی مکانی- زمانی 
بارش پیش‌بینی مدل‌های جوی هستند. مدل‌ها می‌توانند انسجام فیزیکی اجزا مختلف بودجه آب محاسبه شله را 
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حفظ نمایند. همچنین با پیش‌بینی‌های مدل» مطالعه تغییر چرخه آب بر اثر تغییر اقلیم در آینده میسر می‌شود. با وجود 

مزیت‌های فراوان» مدل‌های آب و هوایی در شبیه‌سازی فرآیندهای محلی و بارش مناطق کوهستانی خطا دارند. مقدار 

منطقه‌ای نیز به منطقه مورد مطالعه, تفکیک مکانی و پارامترهای مورد استفاده در آن بستگی دارد (باکس "و برومویچ " 

۶ ). با توجه به مزایا و معایب هر یک از مجموعه داده‌هاء تحقیقات متعددی جهت اعتبارسنجی و مقایسه این 
داده‌ها انجام شده است: 

برخی از مطالعات به مقایسه تصاویر ]11۷1۷ با ساير داده‌های اندازه‌گیری شده پرداخته‌اند (هافمن " و همکاران» 

ماهواره‌ای و سایر داده‌ها مقایسه کرده‌اند. به‌عنوان‌مثال اکتر * و اسلام" (۰۷ ۰ ۷۵ برای پیش‌بینی وضع هوااز حمله 

بارندگی در منطقه بنگلادش از مدل ۷115" استفاده کردند. ایشان این مدل را با دو تفکیک مکانی ۵740 و ۱۵۱۵ 

کیلومترمربع اجرا کردند. هدف از این مطالعه دستیابی به پیش‌بینی بهتر بارش‌های شدید و توانایی هشدار قبل از وقوع 

سیل بود. ات محفقین خروجی مدل 5 را با مجموعه داده 1 ماهواره 1 و باران‌سنج و داده‌های 

رادار مقایسه نموده و به این نتیجه رسیدند که مدل ۷1۷15 با هر دو تفکیک مکانی. از قابلیت خوبی برای پیش‌بینی 

بارش برخوردار است. وونگ" و چیو* (۲۰۰۸: ۱۰۵) به آنالیز ساختار مکانی و زمانی داده‌های اندازه‌گیری شده 

زمینی و اطلاعات استخراج شده از 111101 در هنگ‌کنگ پرداختند. ایشان به این نتیجه رسیدند که تغییرات 

درون‌سالانه زیادی بين اندازه گیری بارند گی و داده‌های سنجش‌ازدور وجود دارد. این پژوهشگران نشان دادند که با 

افزایش اندازه شبکه مقدار میانگین مربعات خطا (18۷1817) و میانگین خحطای مطلق (۷1۸) کاهش می‌يابد. در این 

تحقیق ضریب همبستگی میان بارندگی ۳ ساعته 11011 با اندازه شبکه ۰/۵<۰/۵ و تغییرات درون سللانه زیاد 

۵ به دست آمد. 

احوا"" و همکاران (ع۲۰۱: ۳۷۹ تولیدات 321342 ماهواره ۷۲ را با بارند گی شبه‌سازی شده توسط مدل 

۳ ارزیابی و مقایسه نمودند. منطقه مورد مطالعه ایشان اکوادور و پرو بود. هدف آن‌ها یافتن مناطقی بود که در 

نیت داده‌های مدل و یا سنجش‌ازدور الگوی مکانی - زمانی را بهتر نشان می‌دهند. نتایج حاصل مشسخص کرد که 
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1 در نشان دادن توزیع مکانی میانگین بارندگی سالانه عملکرد بهتری نسبت به مدل دارد. به‌طورکلی هر دو 
مجموعه داده بارندگی‌های ماهانه و سالانه را نسبت به بازه روزانه بهتر تخمین زدند. همچنین در هر دو مجموعه داده 
احتمال تشخیص بارندگی با سرعت‌های کم بیشتر از بارش‌های شدید بود. آزادی و همکاران (۱۳۸۹: 6۱2۸ 
خروجی بارش دو مدل ۷1۷15 و ۷۷18۳ را در مناطق مختلف ايران و یک دوره هشت ماهه مقایسه کردند. نتایج 
نشان دادند که هر دو مدل در مورد بارش فرا پیش‌بینی دارند و هشدارهای نادرست مدل ۷۷۳ بیشتر از 1۷115 
است. هدایتی دزفولی و آزادی (۱۳۸۹: ۱۲۹ پیش‌بینی بارش مدل منطفه‌ای ۷1۷15 را روی ایران مورد بررسی و 
راستی آزمایی قرار دادند. آن‌ها با توجه به ویژگی‌های ناهمواری» اقلیمی و میانگین بارش. کشور را به ٩‏ منطقه تقسیم 
کردند و برای هر منطقه راستی آزمایی جداگانه‌ای را انجام دادند. این پژوهش نشان داد که توانایی مدل برای 
آستانه‌های بارش بسیار کم و زیاد در مناطق جنوبی و شمالی بیشتر است و خروجی مدل در این موارد دفت بیشتری 
دارد. 

پر اساس بررسی منابع مختلف. در ایران تحقیقات اندکی در مورد مقایسه خروجی بارندگی مدل‌های پیش‌بینی 
عددی با محصولات 11311 انجام شده است و اغلب تحقیقات به راستی آزمایی یک مدل یا مقایسه خروجی 
مدل‌های پیش‌بینی عددی با یکدیگر پرداخته‌اند. با توجه به بررسی‌های انجام شده مشخص گردید که اغلب مطالعات 
مذکور از روش‌های آماری برای مقایسات بهره گرفته‌اند. در این تحقیق علاوه بر روش‌های مقایسه آماری» از 
الگوریتم مکانی - زمانی تشابه نیز استفاده شده است. الگوریتم تشابه می‌تواند برای مقایسه کمی همزمان تعداد زیادی 
از نقشه‌ها با چندین متغیر و نیز انجام اين مقایسه به صورت تلفیقی مکانی - زمانی استفاده شود (هرزفلد" و مربام؛ 
۰ ۳). این الگوریتم برای مقایسه کمی همزمان تعداد زیادی داده استفاده می‌شود. چنین داده‌ها می‌توانند در قالب 
نقشه و یا خروجی مدل‌ها باشند. روش مقایسه مکانی با استفاده از الگوریتم تشابه اولین بار برای اکتشاف نفت مورد 
استفاده قرار گرفت (هرزفلد و مریام. ۱۹۹۰ ۶۳). دوربوت " و همکاران (۲۰۰۳: ۱) در بررسی توزیع مکانی تشابه و 
عدم تشابه بین مجموعه داده‌های دما و بارندگی از الگوریتم تشابه استفاده کردند. آن‌ها از این روش برای شناسایی 
ماه‌های مشکل‌ساز در تحلیل داده‌های زمانی استفاده نمودند. ایشان الگوریتم تشابه جبری را به دلیل توانایی‌اش در 
بررسی تعداد زیادی نقشه ورودی انتخاب کردند. 

با توجه به اهمیت داده‌های بارش در پیش‌بینی وقوع مخاطراتی مانند سیل و خشکسالی» همچنین وجود 
محدودیت‌های اقتصادی, جغرافیایی و اجتماعی در افزایش تعداد نقاط اندازه‌گیری اطلاعات بارندگی» یافتن مناطقی 
که در آنجا داده‌های مدل و یا سنجش‌ازدور الگوی مکانی -زمانی بارش را بهتر نشان دهند» ضروری به نظر می‌رسد. 
این امر در مناطق کوهستانی و اقلیم‌های حشک و نیمه‌عشک اهمیت بسزایی دارد. همدف تحقیق حاضر ارزیابی و 
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مقایسه انواع داده‌های بارش با استفاده از روش‌های آماری و الگوریتم تشابه به صورت مکانی و در بازه زمانی 
مختلف به‌منظور استفاده از اطلاعات آن‌ها در نواحی و زمان‌های فاقد داده است. 


۲ منطقه مورد مطالعه 

منطقه مورد مطالعه استان‌های خراسان شمالی. خراسان رضوی و خراسان جنوبی واقع در شمال شرق ايران و بین 
طول ۵۵ تا 7۱ درجه شرقی و عرض ۳۰ تا ۳۸ درجه شمالی می‌باشد (شسکل ۱). اقلیم کلی منطقه نیمه‌عشک تا 
خشک است (احمدیان و همکاران ۱۳۸۰: ۵۹۳). وسعت کل ناحبه ۲۱۳۰۰۰ کیلومترمربع بوده و بلندترین نقطه 


منطقه با ارتفاع ۳۳۰۰ متر بالاتر از سطح دریا در قله بینالود واقع شده و پست‌ترین نقطه با ارتفاع ۲۵۰ متر بالاتر از 


سطح دریا در محل خروجی دشت سرخحس قرار گرفته است. 


۱۱ 
| گراسان شمالی 


شکل ۱. موقعیت منطقه مورد مطالعه 


۳ مواد و روش‌ها 
داده‌های اندازه‌گیری ایستگاه زمینی 
در تحقیق حاضر از اطلاعات ایستگاه‌های هواشناسی سینوپتیکی: ایستگاه‌های اقلیم‌شناسی و ایستگاه‌های آب 
منطقه‌ای خحراسان استفاده شد. اطلاعات مربوط در بازه زمانی ۲۰۰۰ تا ۲۰۱۰ از سازمان هواشناسی کشور و سازمان 
وزارت نیرو دریافت شدند. به‌منظور افزایش دقت محاسبات از میان ایستگاه‌های موجود آن‌هایی که داده‌های گم‌شده 


کی دافعک اتخاب: شنین: 


۱۳۷5 جغرافیا و مخاطرات محیطی شماره بیست و سوم 


داده‌های 0۱۱(۲۸۱۷۲۱۷۲زووز]۱۷ )صمجصمبیعهع]۱۷ 162111 لی1۳0۵1) 


سری هفتم تولیدات 3942 سنجنده 113۷101 به عنوان یکی دیگر از داده‌های ورودی الگوریتم‌های پیشنهادی 
استفاده شد. تولیدات مذکور حاوی داده‌های بارند گی سه ساعته در سطح زمین است و از سال ۱۹۹/۸ در دسترس 
می‌باشند. در این مجموعه داده برای تولید تصاویر بارند گی از ترکیب اطلاعات سنجش‌ازدور (تصویرگر ریزموج, 
رادار بارندگی» پویش گر مادون قرمز و مرئی) و مشاهدات ماهانه محلی استفاده شده است (هافمن " و بولوین " 
کیلومتر) پوشش می‌دهد. در این مطالعه سری زمانی تولیدات 1140۷1۷ از سال ۲۰۰۰ تا ۲۰۱۰ از سایت سازمان 
۸6۵ (/0110://0۳۲-0۵52.20۷) به صورت رایگان دریافت شدند. یکی آر نقاط ضعف داده‌های 711۷۲ 
تفکیک زمانی کم آن است. بطوریکه تعداد عبور آن یک يا دو بار در روز و وابسته به مکان مورد نظر است. این عامل 
سبب از دست رفتن اطلاعات عظیمی در بازه زمانی کوتاه مربوط به طوفان‌ها و رگبارها می‌شود (واند و هیکادان: 
۲۱ ۲ بااین وجود داده‌های 1141۷1 در مناطق کوهستانی که فاقد ایستگاه‌های باران‌سنجی می‌باشند. مهم 
بوده و می‌تواند در مطالعات هیدرولوژیکی و هواشناسی مورد استفاده قرار گیرند. 
داده‌های پیش‌بینی مدل 1۷۲۷5 
با پیشرفت علوم کامپیوتر و هواشناسی سینوپتیکی. شمار زیادی از مدل‌های عددی میان‌مقیاس برای شبیه‌سازی و 
پیش‌بینی حرکات جو و متغیرهای هواشناسی استفاده شدند. مدل ۷۲۷۲5[ (نسل پنجم مدل میان‌مفایس ۱۸۹/۲۵۱۸ 
عاقملا) نیز از این قاعده مستننی نبوده و برای پیش‌بینی متغیرهای جوی در مقیاس سینوپتیک به کار می‌رود. این مدل 
از معادلات اولبه غیرهیدرواستاتیک برای پیش‌بینی‌ها استفاده می‌کند. محتصات عمودی مدل 0 و مخحتصات اققی آن 
4تع- ۸۵۵۷۵ می‌باشد. ب‌طور حلاصه مدل 1۷1۷/5 یک مدل منطقه محدود. غیر ایستا و متأثر از عوارض زمین 
است که برای شیه‌ساژی با بستن‌بیتی گردشن میان مقیاس جو طراحی شده است. همچنین شبکه‌بندی مدل 5 بر 
اساس کاربردهای خاص قابل تنظیم انتیت (رلا و همکاران» ۵ در ان تحفیق داده‌های مربوط به پیش‌بینی 
مدل 1۷15 از سسابت سس از مان 130۸۸ ( 01105://07۷۷۷۷۰0۵06۰1022.20۷) و سس زمان ۲۲۸۶ 
(۳1105://۷۷۷۷:/۰0627.60/) دریافت شدند. اطلاعات دریافت شده تصاویری با تفکیک مکانی ۰۱/۵۰/۵ درجه و دوره 
زمانی ۲۰۰۰ - ۲۱۱۰ بودند که بر اساس طول و عرض منطقه مورد مطالعه تهیه شدند. به دلیل وجود خحطاهای 
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سال ششم ارزیابی داده‌های بارند گی حاصل از ماهواره 15۷1۷ ... ۱۹۹ 


قاعده‌مند در برونداد مدل‌های پیش‌بینی عددی. قبل از استفاده از داده‌ها به عنوان ورودی الگوریتم‌ها؛ پس‌پردازش 
آماری روی خروجی مدل انجام شد. علاوه بر اطلاعات ذکر شده داده‌های مربوط به توپوگرافی از سازمان 
نقشه‌برداری زمین‌شناسی ایالات متحده (7505) (110://120060۷76.0585.207) با تفکیک مکانی ۳۰ انیه دريافت 
شبکه‌بندی داده‌های زمینی بارندگی 
داده‌های ماهواره و مدل‌های حوی به صورت تصاویر شبکه‌بندی شده ارائه می‌شوند در حالی که مشاهدات 
زمینی اطلاعات نقطه‌ای فراهم می‌آورند. یکی از مشکلات موجود در زمینه مقایسه داده‌های به دست‌آمده از منابع 
مختلف» نمونه‌برداری مکانی متفاوت آنهاست؛ بنابراین داده‌های نقطه‌ای باید با استفاده از روش‌های درون‌یابی مکانی 
به صورت تصاویر شبکه‌بندی با تفکیک مکانی یکسان با ساير داده‌ها تبدیل شود. در مورد انتخاب مناسب‌ترین روش 
درون‌یابی مخحصوصاً در مناطقی با توپوگرافی پیچیده مطالعات فراوانی صورت گرفته است. از میان آن‌ها روش 
زمین‌آمار به عنوان روشی دقیق و قوی در این زمینه شناخته شده است رکوازت » 4 مطالعات بسیاری اثبات 
کرده‌اند که برای بارندگی روش کربجینگ با روند خارجی («اتل16) که از ارتفاع به عنوان متغیر ثانویه استفاده می‌کند 
در مناطق با توپوگرافی پیچبده بهتر از سایر روش‌ها عمل می‌کند (میسون " و فری * ۲۰۱6 1۵8۳؛ توبین * و همکاران, 
۱ با وجود اين. در این تحقیق الگوریتم‌های دیگر زمین‌آمار نیز بر روی داده‌ها آزمایش شد و براساس نتایج 
اعتبارسنجی (یعنی مقادیر میانگین مربعات خطا (۷1517) که با استفاده از معادله ۱ محاسبه شدند) از روش 16۳۲ 
برای درون‌یابی بارند گی استفاده شد. معادلات مربوط به روش‌های زمینآمار توسط وبستر" و اولیور*(۲۰۰۱) به 
تفصیل شرح داده شده‌اند. الگوریتم زمین‌آمار با استفاده از ترم افزارهای ب+وی) و 20( بر روی داده‌ها اجرا شدند. 
پس از تهیه نقشه‌های نوزیع مکانی بارندگی» از دو روش جهت مقایسه داده‌ها استفاده سل روش اول استفاده از 
تصاویر بارندگی ماهواره 1191۷1 ريشه میانگین مربعات خطا (۷]5۳() در مقایسه با مقادیر داده‌های زمینی 
بارند گی برای دو محصول ذکر شده. محاسبه شد (معادله 0 روش دوم استفاده از الگوریتم تشابه و جستجوی مکانی 
تشابهات و تناقضات میان داده‌ها بود. معادلات مربوط به الگوریتم تشابه در ادامه توضیح داده شده‌اند. 
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در این روابط ,2 مقادیر پیش‌بینی شده بارندگی در نقطه و ز2 مقدار اندازه‌گیری شده بارندگی است. 1 تعداد 
کل داده‌ها می‌باشد. 

الگوریتم تشابه 

به‌منظور یافتن مکان‌هایی که در انشا مدل ۷/5( و داده‌های سنجش‌ازدور الگوی مکانی- زمانی بارش را بهتر 
نشان داه‌اند از روش الگوریتم تشابه استفاده شد. در این الگوریتم مقادیر تشابه از معادله (۳) به دست می‌آیند. اگر 1 
نشان دهنده تاحبه نقشه و 1 تعداد نقشه‌های ورودی و () ,1۳ مقدار استاندارد شده نقشه 3 در مکان باشد» 
برای هر نقطه شبکه ( 26/۷ المان‌های ماتریس اختلاف 0(67۳)( با استفاده از معادله (۵) شکل می‌گیرند(مریام! 


و جوت ۱۹۸۸: :)٩‏ 
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در این روابط 6 تعداد حالت‌های مقایسه دو به دو در میان 7 نقشه ()ءو2 المان‌های ماتریس اخحتلاف 
۴( برای هر نقطه شبکه ( 561 ()و7 مقدار استاندارد شده نقشه و11 در مکان :۰2 (2)ع70 مقدار 
اتاندارد شده نقشه 1/2 در مکان ۰26 ۰,1 1,۰ < 5و 9,. ..,1 < ]می‌باشند. با حل ماتریس اختلاف برای هر 
نقطه شبکه. میزان تشابه و تناقض میان تصاویر آشکار می‌گردد. (0) 7 میانگین اخحتلاف بین نقشه‌ها در مکان 4 را 
نشان می‌دهد و بنابراین اگر نقشه‌ها مشابه باشند مقدار (2) ۳ کم است و اگر متفاوت باشند اين مقدار زیاد است. در 


این تحقیق قبل از اجرای الگوریتم تشابه با استفاده از رابطه زیر استاندارد سازی انجام گرفت: 
0 ز۷771۳0 - ززلا 


ز2 - 177102 < (رلا)مت 
در رابطه )1( زا7171 و 1717102 مقادیر بيشینه و کمینه متغیر مشاهده شده را می‌باشند. بطوریکه < ز 


ز۰..,1 ,1 تغییر می‌کند و 7 تعداد مشاهدات متغیر ( است (هامان " و هرزفلد* ۱۹4۱ 1۲۱). در این پژوهش از 
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تبدیل خطی محدوده داده‌ها با فواصل (۰,۱) استفاده شد و سپس محاسبات در داخل منطقه مورد مطالعه انجام 
گردید. پس از عملیات استانداردسازی مقادیر (6) ] برای تمام منطقه مورد مطالعه محاسبه و در هر موقعیت 22 
یک مدل شبکه‌ای برای نقشه تشابه ۳ ساخته شد. همچنین در مورد تحلیل کیفی و کمی تشابهات و تناقضات موجود 
در نقشه‌های خروجی الگوریتم تشاب به صورت زیر عمل شد: 

چنانچه نسبت محدوده استانداردسازی به کار رود آنگاه مقادیر ۳ در نقشه تشابه در بازه صفر و یک قرار 
می‌گيرند. در این صورت مقادیر نزدیک صفر نشان دهنده تشابه خوب بین نقشه‌های ورودی در یک مکان خحاص 
بوده و مقادیر بالات نشان‌دهنده تشابه ضعیف‌تر می‌باشند. با توجه به مطالعات انجام شده در مقایسه ٩‏ نقشه یا بیشش 
اگر مقدار ۳ زیر ۰/۲ باشد تشابه خوب و بالای ۰/۵ تشابه ضعیف را نشان می‌دهد (مریام و جووت. 1۹۸۸). 
طبقه‌بندی مقادیر تشابه (۳) در جدول (۱) نشان داده شده است. در این مطالعه الگو ریتم‌های استانداردسازی و 
تحلیل‌ها با برنامه‌نویسی در نرم‌افزار ۷۸۲1۸3 برای هر نقطه شبکه از نقه‌ها اجرا شدند. برنامه الوریتم تشابه نیز 
در نرم افزار ۷1۵1۸ نوشته و بر روی تمامی نقشه‌ها اجرا شد تا تشابهات بین مجموعه داده‌ها و نتایج مدل از 
لحاظ تحلیل مکانی و زمانی کشف شود. در نهایت توسط نرم افزار ]۷۸ نتایج حاصل ترسیم شدند. 


حدول ۱. طبقه‌بندی مقدار ۲ 


میزان تشابه م۷2 ۲ 
خوب ۰/-۱/۱۱ 
اشنا توت ۱/-۰/۲ 
متوسط ۱/۲۲ 
تا میت ۰/۳ 
ضعیف ۰/۵ 


۶. بحث و نتایج 

زمین آمار 

تلو ۲ افارشش وهای اه یدای آمان زا رای وهای ۱ فان مر ده 
مقافزر دا ان مس هید که با سمل زوش عا ان اند کی رای ک‌داتن مره بعتار اقااز یه 
صفر نزدیک‌تر باشد» دقت روش در پیش‌بینی مقادیر بیشتر است. طبق جدول (۲) در شبیه‌سازی مقادیر بارندگی؛ 
تقاوات ین خفت دون قر قضی ار مایها مس داز تست اما تطی ان زو اتعاوی توا کشری تس 
به روش 001 دارد. در اين مطالعه با توجه به مقادیر ماهانه 18 روش 16710 با متغیر ثانویه مدل رقومی ارتفاع 
منطقه (0۳۷1) برای درون‌یابی بارندگی و تحلیل مکانی بارندگی انتخاب شد. 


۱۷۳ جغرافیا و مخاطرات محیطی شماره بیست و سوم 


جدول ۲. مقادیر 115۳ روش‌های 15 و >01) برای تهیه نقشه‌های توزیع مکانی بارندگی در مقیاس ماهانه و 
سالائه 

امه | 0۵ | ۱0۵۲ | 060۲ ز ظاک | متام | تال | ها | تخشالا [ ۸۳۴ ۱ ۲۸ | 1۳۳۵ | از | بارندگی 

۵٩‏ ۱ | ۲ | ۰۸۱۵ ]۱۸۲ ۰۲ ]۸۵| ۰۲۵ ۱۵ ۱۲۱۲ ۰۱۴ ۰۲۵ ول 


01 ۱/۳ ۱۰/۶ | ۰/۶۷ | ۰/۳۷ | ۰/۵۳ ۱۰۸۷۸ ۱ ۰/۱۷ | ۰/۲۳ ۱۰/۱۷ | ۰/۳۶ ۲ ۲ ۱ ۰/۷۷ 


اعتبارسنحی داده‌های ۷۲1۷15 110/۷۷" و داده‌های شبکه‌بندی شده ایستگاهی در مقیاس ماهانه 

شکل (۲) مقادیر 121۷161 مدل 1۷1۷15 و 11۷1*1 را در مقیاس زمانی ماهانه نشان می‌دهد. مقادیر 115۳ با 
توجه به میانگین ماهانه داده‌های بارندگی ایستگاهی به دست آمده‌اند. هر چه مقدار 11۷15۳ کمتر و به صفر 
نزدیک‌تر باشد تخمین‌ها بهتر انجام گرفته‌اند. با توجه به شکل (۲)» مدل ۷1۷5 در ماه‌های گرم مقادیر 1315۳۴ 
بیشتری دارد. بدین معنی که بارندگی در ماه‌های گرم با دقت کمتری نسبت به ماه‌های سرد پیش‌بینی شده است. 
همچنین بیشترین مقدار مربوط به اواخر فصل بهار و اوایل تابستان می‌باشد. دوربوت " و همکاران (۲۰۰7: ۱) نیز در 
بررسی‌های خود به نتیجه مشابه دست یافتند و دلیل این امر را بارندگی‌های همرفتی و عدم توانایی مدل در پیش‌بینی 
دقیق‌تر اين نوع از بارندگی عنوان کردند. الگوی تغیبر مقادیر 18۷1517 در نمودار شکل (۲) برای 1180۷1۷1 معکوس 
مدل 115 است. بیشترین مقدار 13۷157 برای 118101 در ماه ژانویه و تقریباً برابر با ۱۶ به دست آمده است. دلیل 
این امر ممکن است مشابهت پراکندگی انرژی مایکروویو رفتار بخ روی زمین مانند بارش و یا باران یخ زده در 
اتمسفر باشد. نتایج اچوا" و همکاران (۲۰۱6 ۳۱۷۹) نیز مانند یافته‌های تحقیق حاضر نشان دادند که تولیدات 
ماهواره‌ای 1180۷1۷1 در مناطق پوشیده از بخ و ماه‌های سرد فصل زمستان عملکرد ضعیفی دارند. در ماه‌های گرم؛ 
سنجنده ]1181۷ بارندگی را بهتر از مدل برآورد کرده بود؛ اما با توجه به محدوده تغییرات 120۷1۹۳ مدل ۷15[ 
مقادیر بارندگی را با نوسان ماهانه کمتری نسبت به 18۷11 پیش‌بینی کرده است. به جز ماه‌های تابستان در بقیه 
سال. مقادیر :16۷15 مدل کمتر از 116۷101 می‌باشد. بدین معنی که مدل ۷1۷15 بارندگی را با دقت بیشتری نسبت به 


1 پیش‌بینی کرده است. نتیجه به دست آمده مشابه نتایج اکتر" و اسلام * (۲۰۰۷: ۷۵) بود. 
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زمان 


شکل ۲. اعتبارسنجی داده‌های ۱1315 و 118۱131 در مقیاس ماهانه 


مقایسه مکانی-زمانی داده‌های اندازه‌گیری ایستگاه زمینی. 1۷115 و 116۷1۷1 با استفاده از الگوریتم تشابه 

در بکارگیری نتایج الگوریتم تشابه دو عامل یعنی مقدار تشابه (۳) به دست آمده از معادله (۳) و یکنواختی توزیع 
مکانی مقادیر تشابه در منطقه مد نظر قرار می‌گیرد. بدین صورت که هر چه مقعدار ۳ کمتر و توزیع آن یکنواخت‌تر 
باشده تشابه مکانی» بیشتر است. نوع تشابه در هر ناحیه نیز از جدول (۱) تعیین می‌گردد. 

تغییرات درون سالانه داده‌های مشاهده‌ای بارندگی (نقشه‌های حاصل از 16۳۴1 

شکل (۲) تغییرات درون سالانه بارندگی داده‌های مشاهده‌ای را نشان می‌دهد. شکل ۳-لف مربوط به فصل سرد 
و شکل ۲سب مربوط به فصل گرم می‌باشد. با توجه به شکل مشاهده می‌شود که در فصل سرد مقادیر ۳ در منطقه 
مورد مطالعه بیشتر بین صفر تا ۰/۱ تغییر داشته و در سطح وسیعی از منطقه مقادیر تشابه بین ۰/۱-۰/۰۵ می‌باشند. 
همچنین مقادیر تشابه در ناحیه مورد مطالعه با توزیع یکنواختی پراکنده شده‌اند که این نشان دهنده تغییرات درون 
سالانه کم بارندگی در فصل‌های سرد سال است (شکل ۳- الف). در شکل (۳- ب) که مقادیر بارندگی برای 
فصل‌های گرم سال می‌باشد. مقادیر تشابه در جنوب و شمال منطقه بین ۰/۱ تا ۰/۲ تغییر می‌کند که نسبت به همین 
مقادیر در فصل‌های سرد سال بیشتر است و داده‌ها تغییرات درون سالانه بیشتری نسبت به فصل سرد دارند. همچنین 
نواحی مرکزی منطقه مورد مطالعه از مقادیر تشابه بالاتری برخوردار بوده و مقادیر ] در این ناحیه بین ۰/۲ ۰/۳۵۱ 
تغییر کرده است و بیانگر این است که تغییرات درون سالانه بارندگی در فصل‌های گرم سال و در نواحی مرکزی 
بیشتر بوده است. یکی از دلایل این تفاوت می‌تواند وجود خطا در داده‌های مشاهده‌ای باشد (دوربوت " و همکاران 
۰۹ با توجه به شکل ۳ - الف و ۳ -ب می‌توان گفت که تغییرات درون سالانه بارندگی در فصل گرم خیلی 
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بیشتر از فصل سرد بوده است. تغیبر نوع بارندگی و همرفتی بودن آن می‌تواند دلیل این امر باشد. به خصوص در 


نواحی مرکزی که دارای ارتفاع بلندتر و توپوگرافی پیچیده‌تری است. 


شکل ۳. تغییرات درون سالانه داده‌های مشاهده‌ای برای متغیر بارندگی. الف : فصل سرد. ب : فصل گرم 


تغییرات درون سالانه خروحی مدل ۷۲5 

تغییرات درون سالانه بارندگی مدل 1۷1۷5 در شکل (4) نشان داده شده است. با توجه به شکل, این تغییرات در 
مدل ۷1۷85 کم بوده و در هر دو نقشه تشاب فصل سرد (شکل 4- الف) و فصل گرم (شکل ۶-ب) با مقادیر ۳ 
قابل مشاهده می‌باشند. با توجه به شکل مشاهده می‌شود که در فصل سرد مقادیر "1 در منطقه مورد مطالعه بیشتر بین 
۵ ۷ تغییر داشته و اغلب مساحت منطقه را مقادیر تشابه بین ۰/۱۵-۰/۱ در برگرفته است. همچنین 
پراکندگی مقادیر تشابه در منطقه از توزیع یکنواختی برخوردار می‌باشد. علاوه بر این مشاهده می‌شود که نقشه‌های 
تشابه مدل 1۷1*15 بین تابستان و زمستان متفاوت است. با توجه به شکل ۶- ب نواحی مرکزی منطقه مورد مطالعه از 
مقادیر تشابه بالاتری برخوردار بوده و مقادیر ۲ در این ناحیه بین ۰/۱۰/۱۵ تغییر کرده‌اند. این موضوع نشان 
می‌دهد که تغییرات درون سالانه بارندگی در فصل‌های گرم سال و در نواحی مرکزی بیشتر بوده است. ممکن است 
وجود ارتفاعات و توپوگرافی پیچیده این ناحیه دلیل آن باشد. با توجه به جدول (۱) که طبقه‌بندی مقادیر ۳ را نشان 
می‌دهد و با توجه به مقدار "1 به دست آمده در شکل. به‌طورکلی مدل بارندگی فصل سرد را حوب و بارندگی فصل 
گرم در نواحی مرکزی منطقه را نسبتاً ضعیف پیش‌بینی کرده است. نتیجه اخبر با تحقیق انجام شده توسط هدایتی 


1 نك ارزیابی داده‌های بارند گی حاصل از ماهواره ۷ ... ۱۷6۵ 


دزفولی و آزادی (۱۳۸۹: ۱۲۹) مطابقت داشت. ایشان چیرگی نقش نایایداری کزفشاری نسبت به ناپایداری همرفتی 


در سامانه‌های فصل زمستان را دلیل موفقیت مدل در پیش‌بینی‌های فصل سرد عنوان کر ده‌اند. 


۳ ۲ 


شکل . تغییرات درون سالانه خروحی مدل 5 برای متغیر بارندگی. الف : فصل سرد. ب : فصل گرم 


تغییرات درون سالانه بارندگی ماهواره 11۷11" 

شکل ۵ الف و ب میانگین مقادیر تشابه محاسبه شده در سری زمانی ۲۰۱۰-۲۰۰۰ داده‌های ]1181۷ را نشان 
می‌دهند. با توجه به شکل ۵- الف مشاهده می‌شود که در فصل سرد. جنوب منطقه دارای مقادیر تشابه ۰/۰۵ -۰/۱ 
است و توزیع مقادیر نیز یکنواخت می‌باشد؛ اما نواحی شمالی و مرکزی, از مقادیر تشابه بالاتری برخحوردار هستند. 
بطوریکه مقدار ۳ در این نواحی از ۰/۲ تا ۰/۲۵ متغیر است. همچنین توزیع مقادیر چندان یکنواخت نیست. با توجه 
به این نکات می‌توان نتیجه گرفت که تولیدات 118۷1۷1 مقادیر بارندگی در فصل سرد را در جنوب منطقه نسبتا 
خوب و در شمال و مرکز آن نسبتاً ضعیف برآورد کرده‌اند. تغییرات درون سالائه فصل گرم نیز در شکل ۵-ب نشان 
داده شده است. چنانچه از شکل مشخص است بیشتر مساحت منطقه را مقادیر ۰/۱ تا ۰/۲ در برگرفته‌اند و توزیع 
مقادیر تشابه در منطقه یکنواخت‌تر از فصل سرد (شکل ۵-لف) می‌باشد. بیشترین مقدار ۳ در قسمت شمال شرقی 
منطقه مورد مطالعه مشاهده می‌شود؛ اما به‌طورکلی 118۷1۷1 بارندگی فصل گرم را نسبتاً حوب برآورد کرده است. 
همچنین با مقایسه شکل ۵-لف با شکل‌های 6-الف و ۳-لف مشخص می‌شود که رفتار مقادیر بارندگی ]7101۳ 
عکس داده‌های مشاهداتی زمینی و 1۷1۷15 است؛ مانند مطالعه وونگ " و چیو" (۲۰۰۸: ۱۰۵) نتایج این تحقیق نیز 


۹۵ 1 
الط 2 
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نشان دادند که به‌طورکلی تغییرات درون‌سالانه زیادی بين اندازه گیری بارند گی و داده‌های سنجش‌ازدور وجود دارد. 
همچنین با مقایسه شکل‌های ۳ ۶ و ۵ با یکدیگر متوجه می‌شویم که مدل 5 بهتر از 11*۷1 تغییرات فصلی 


بارندگی را نمایان کرده است. چرا که تغییرات آن در هر دو فصل, با داده‌های مشاهدات زمینی مشابه است. 


شکل ۵. تغییرات درون سالانه تصاویر بارندگی 71 الف : فصل سرد. ب : فصل گرم 


۵ نتیجه گیری 

در این مطالعه داده‌های بارندگی مدل 4۷115 سنجنده 100۷11 و مشاهدات ایستگاهی در محدوده سال‌های 
۰۱۰-۰ و شمال شرق ایران توسط روش آماری و الگوریتم تشابه مقایسه و صحت‌سنجی شدند. هدف از این 
کار بررسی امکان استفاده از داده‌های مدل‌سازی شده در مناطق فاقد داده و شناسایی مکان‌ها و زمان‌هایی بود که 
محصولات مذکور با یکدیگر و داده‌های واقعی همخوانی بیشتری دارند. نتایج اعتبارسنجی آماری به صورت ماهانه 
نشان دادند که مدل 1۷115 در ماه‌های سرد بهتر از ماه‌های گرم مقادیر بارندگی را پیش‌بینی کرده است. نتیجه به 
فقس آمده با نتایج دوربوت " و همکاران (۲۰۰7: ۱) تطابق داشت؛ اما تصاویر سنجش‌ازدور بر عکس مدل ۷1۲5 در 
ماه‌های سرد بیشترین خطا را داشتند. با توجه به این مورد توصیه می‌شود کاربرد داده‌های 116۷1۷1 در فصول سرد و 
نواحی پوشیده از برف و یخ با احتباط انجام شود. همچنین در اغلب موارد مدل 1۷1۷15 در پیش‌بینی بارندگی 


اون داشت. آزادی و همکاران (۱۳۸۹: ۱۶۸) نیز در بررسی‌های خود به نتیجه مشابه دست يافته‌اند. بر خلاف 
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یافته‌های اچوا! و همکاران (۲۰۱6: ۳۱۷۹) که در مناطق اکوادور و پرو انجام شد. در اين مطالعه (شسمال شرق ایسران»؛ 
داده‌های 1180۷111 نسبت به خروجی مدل پیش‌بینی عددی دارای حطای بیشتری در برآورد بارندگی بود. 

در این پژوهش مجموعه داده‌های مختلف تغییرات بارندگی را در منطقه مورد مطالعه به‌حوبی تخمین زده بودند؛ 
اما برآوردهای مناطق مرکزی مخصوصاً ارتفاع بالای ۱۵۰۰ متر در محصولات دقت پایینی داشتند. در بررسی تغیسرات 
درون سالانه و فصلی. مشخص شد که تفاوت زیادی میان تغییرات درون‌سالانه تصاویر سنجش زدور و مقادیر واقعی 
بارندگی (1۳۳) وجود دارد؛ اما تغییرات درون سالانه مدل با داده‌های واقعی همخوانی داشت. 

نتایج تحقیق حاضر نشان دادند که در صورت نبود اطلاعات ایستگاهی کافی برای بارندگی در یک منطقه 
می‌توان از داده‌های مدل ۷1۷5 و 112۷11 برای پر کردن خلا موجود در این مورد استفاده کرد. این در حالیست که 
پیوستگی مکانی-زمانی داده‌های بارندگی در مدل‌سازی هیدرولوژیکی و مدیریت هواشناسی کشاورزی از اهمیت 
بالایی برخوردار است. در موارد کاربردی مانند هشدار سیل یا مدیریت مخاطراتی مانند خشکسالی» نیاز به داده‌هایی 
با تفکیک زمانی و مکانی بالا وجود دارد. با توجه به کمبود امکانات و نبود ایستگاه‌های باران سنجی و هیدرولوژی؛ 
استفاده از داده‌های مدل‌سازی بارندگی و محصولات سنجش ازدور می‌تواند راهگشا باشد؛ اما یافته‌های تحقیتی حاضر 
نشان دادند که نمی‌توان فقط از یک مجموعه داده برای تمام مناطق و یا فصول مختلف در نواحی فاقد داده استفاده 
کرد؛ بنابراین پيشنهاد می‌شود در آینده علاوه بر صحت سنجی مکانی- زمانی به جای روش‌های ساده آماری» 


روش‌های تلفیق داده نیز برای تولید داده‌های با کیفیت و دقت بالا مورد بررسی قرار گيرند. 
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زمین و فضا. ۳۱ (۳): ۱۱۵ تا ۰۱۲۹ 
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.1295۰ 

01 ۱۲21۴ 00- 1ج طا ممتاهاامتهعزم ۵ فلورولقصه لهتاهم5 .2014 وین بلع۳1 عک .ثم ,1۷128501 

۰ ۱۷۱۵۵۵ ممامتاهتاه مه فماممم منصمهتهممه ملهعو-تاانا صز فلمصاعصظ ممتعمآمرم 

45343۰ ,18 ,90 9۷5۲ طات۲ 

ما باتاطه۱۷]6 رل ب۷1۵065۵ ب.ظ ماعطع6992ظ و.ظ مفتصتاط 1 بن) وتع۷0۷۵11) ببظ ب02۵[12 ,.بابه ب,زاع]1۷ 

مه مها امن و0۲ رع9۵)60۲0[0 ۱۷۷ لهج ۷۱۷۲5 ۵۶ )مهمرحظز ۰1۳6 .2009 ,۲.۷ رتع6۱0کز 

6611-3۰ ,(9)17 ,یا ۵ اه و۸2۳۵ .عصمتاهلنهاف ام0مهه 0۳۱۲۱۴۱۴۲ 

5 ۱/۲۲۵۲۶ 00۳9۵211808۰ رفح متقصصفطا ۵۶ فلمطام۷ .1988 ,.ت) .۲۰ ,169/۵0 عک .۲ .1 مطه۱۷۱۵۳ 
9-۰ ,۳۲655 حصتاصماظ ویلظ محطه۱۷]۵ ۳۰ 1۳۰ رقع0760۱۱۵1۷6۳۱۱۵1) 77 

و.لا] ۷۰ بامففوتارظ بعنا .لا بهمصمنان ویو رمطه]2۷۵) ورظ مقتطم۱1 ۱۷۷ ویظ ,0و۵ وم م602ظ۱ظ بنظ م0602 
۳ باحصا ممتامالنآمع0 ۰ 3۳842 ۰۲۳۴۵۷۲۸۷۲ 0 ممتهنطله7 . ,2014 
سچوما۵۲۵ ۳۵۲۱۵۰ 20 ۵0۲۵00۲ظ ۵۶ مملوع صحعع0مه - ۳۵۵1110 معط تمعن ممتاهانصته ممتاماام2۳60 
,3179-3۰ و6 ,۹۵6۱668 5۵۹/۵۱ ۲۵۲1 0۳0 

و.گ ۲۰) بطفحطن ۲۲ 20 ر.ط و1۱۷6 نا رتقلل۷1 وماصهو ری بامععنا۲۳ ,۷ رته‌تطامط ۷ .با ۷ بلهعظ90 
مد ممصمصمعن (۸ظ(1) ولو راممه ممتاماانز۳:۵۵ م)۱۵۱۱۵۵۵/۵[11 ۲۳۴۱۸۷۲ ۵۶ ومتلهتلم۲ظ .2011 
۰ رما نامع لهز0مصعا له لفتاهمود وه «ممع0معمع فلا هه م210ع دص تقتامعی عمط 
2049۰ ,13 ,90 ,)9۷۹ طاتو۲ 

فص صصهتعم۲ 00۸ )مومع ممتاهازن۵۵:ظ تحطامات) تقممتمصمهماه۱ .2007 ریق .ظ رطاتظو 
۲۳۳۵0 مصتنناوهع]۱۷ :.۲ .۲ پلتتا 1 هه بیط رتفتاحظ و۷ متطه1۵۷1272 1 ,مهن صظ تصمتفو1 
-611 ,28 ,۴۵:62 ععمفط 0۱۵92 مد قعمصهب 4 .ععتاباظ معط 20 ۱۲9۸1] نا :عمهمو م۶۲ 
0۰ ,120019 مها 1۳0۲0۲60 رتم‌ومتنجد بددم 

0 .2011 یه م,100فص 24 بیظ مع0تظ و.ظ ۷ مهعصقاته۲ ببا مقصتاهع1 بن) مطلفا 1[ 
۵ ۳/19 2 طا مممتاهااممه معمام بط تم مهم لهماهما0۲0عا۵ ۵۶ ممته1۳۵1۵01 
77-۰ ,401 ,۳۱۷0۲01 .1 .۲۵۵100 


سال ششم ارزیابی داده‌های بارند گی حاصل از ماهواره 15۷1۷ ... ۱۷۹ 


5 0۲۵۵101210 ۵ مملاهن2 ۳۷ :2011 یط رتعاه۷۷۵ 200 وبا بله62۷6 ,۷۷.۰ با۱]261ظ ر.ظ ,27 ۷۷ 
,16۵50۲۰ 2۸۵۲ ۷۷ .عم ۵9۵۵0۷۵ 16501۲۵۵5 ۷۷/۵۵۲ 2 بصلهه) فتامصه) رامق مامرممن 0۷۵۲ 
1222-10۰ ,34 

مطامل :ه)ععطامنطن) .مافتامهمتهه له)مممصمم وه ۲0۲ وما060۹]209 ,2001 وی ,۲ و01۷۵ عک وی ,۱۷۷60۹61 
0۰ 6 ۷۷116۷ 

0ص مادنا معباهت) متفگ ۵۶ فلور امه ۱68۵0۲2 لمح ۵22 .20086 ویو بسا متاللن ‏ رل .۲ ۷۷۰ ,۷۷۵82 
-105 و(4)2 رععع0۵ ۱0۳۱۱۱۵۵۸ 0608۳۵۵6 5۵8۵۰ ۳۱۵9۵ 1 ولج۷عنام [املمتی۳۴ 1۳۴/۷۲۷۲ 
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زار ر‌ محیطی 


شماره بی 


